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SECTION  1 


INTRODUCTION  AN!)  SUMMARY 

The  objectives  of  this  program  are  to  study  the  physics  of  microwave 
and  millimeter-wave  semiconductor  devices  under  optical  illumination,  to 
study  the  techniques  of  modulating  semiconductor  lasers  at  microwave 
frequencies,  and  to  combine  the  above  two  aspects  to  achieve  optical 
control  of  microwave  and  millimeter-wave  semiconductor  devices.  During 
the  second  quarter,  we  completed  calculations  on  the  feasibility  of 
optical  injection  locking  of  millimeter-wave  1MPATT  oscillators,  con- 
tinued with  the  optical  injection  locking  experiment  of  X-band  Si  1MPATT 
oscillators,  and  carried  out  calculations  on  the  mode-locking  charac- 
teristics of  external-cavity  GaAlAs  injection  lasers. 

The  feasibility  of  achieving  optical  injection  locking  of  millimeter 
wave  IMPATT  oscillators  was  investigated  by  calculating  the  efficiency 
of  subharmonic  optical  injection  locking  of  IMPATT  oscillators.  The 
calculations  were  carried  out  by  adding  a time-varying  reverse  satura- 
tion current  term,  which  was  generated  by  the  incident  optical  signal, 
to  the  IMPATT  avalanche  equation.  Solving  the  equation  yielded  a diode 
external  current  component  related  to  the  injected  signal.  This  current 
was  then  taken  as  the  injection  source  in  the  oscillator  equivalent  cir- 
cuit to  calculate  the  locking  range.  The  results  indicated  that  a 
locking  range  on  the  order  of  100  MHz  could  be  achieved  for  X-band 
IMPATT  oscillators  if  modulated  optical  signal  was  converted  into  rf 
photocurrent  efficiently. 

The  injection-locking  characteristics  of  X-band  Si  IMPATT  oscil- 
lators were  further  investigated  by  measuring  the  locking  range  as  a 
function  of  locking  gain.  Subharmonic  locking  was  used  with  a frequency 
ratio  ranging  from  3:1  to  8:1.  A locking  range  of  several  megahertz 
was  achievable  with  an  oscillator  at  8.11  GHz.  Even  subharmonic  locking 
was  preferable  to  odd  subharmonic  locking  for  our  oscillators. 

The  transmission  characteristics  of  two  different  Fabry-Perot 
resonators  were  studied  to  find  an  optimum  external-cavity  configuration 
for  the  injection-laser  mode-locking  experiment.  The  first  resonator 
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1 


consisted  of  an  injection  laser  of  length  P.^  and  a mirror  located  at  a 

distance  P.  from  the  laser.  The  inner  facet  of  the  laser  was  coated 
o 

with  an  antireflection  (AK)  coating.  Therefore,  this  combined  resonator 

acted  as  a simple  Fabry-Perot  resonator  with  modes  that  are  regularly 

spaced  in  frequency  by  Af  = c/2(P,  + nP.,),  where  n is  the  index  of 

o 1 

refraction  of  the  laser  medium.  The  second  resonator  also  consisted  of 

an  injection  laser  and  a mirror;  the  laser  facet,  however,  was  not 

coated.  In  this  case,  we  actually  had  three  sets  of  Fabry-Perot 

resonators  with  lengths  P.  , V,,  + Si  , and  P.,  . Therefore,  the  mode  struc- 

o 1 o 1 

ture  was  more  complicated,  consisting  of  three  sets  of  modes  with  fre- 
quency spacing  c/2Z^,  c/2  (P^  + nP,^),  and  c/2nP,^,  respectively.  However, 
since  ii.  >>  P,^,  the  first  two  sets  of  modes  would  have  roughly  the  same 
spacing  and  be  located  alongside  one  another.  Since  the  second  reson- 
ator was  found  to  have  a higher  Q factor  than  the  first,  it  might  be 
easier  to  achieve  mode-locking  in  this  resonator. 
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SECTION  2 

CALCULATIONS  ON  OPTICAL  INJECTION  LOCKING  OF  IMPATT 

OSCILLATORS 


A.  BASIC  EQUATIONS 

The  equation  that  describes  the  avalanche  current  density  J of  an 
IMPATT  diode  is  given  by* 

+ Js  , (1) 


a dJ 
2 dt 


li 

I 


a dx  - 1 


where 


T 

a 


£ /v  is  the  avalanche  zone  transit  time 
a 

£ is  the  avalanche  zone  length 
a 

(X  is  the  impact  ionization  coefficient 

J is  the  reverse  saturation  current  density 
s J 

v is  the  saturated  carrier  velocity. 


The  ionization  coefficient  a is  a function  of  the  electric  field  in  the 
avalanche  region  and  can  be  approximated  by  the  following  expression: 


u = 


(m  = 6 for  silicon) 


Assume  that  the  electric  field  distribution  inside 
as  shown  in  Figure  1,  where  the  field  varies  from 
the  drift  region  value,  E^ , with  dE/dx  constant, 
that 


the  IMPATT  diode  is 


its  peak,  E^,  down  to 
Therefore,  it  follows 
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0 


AVALANCHE 

REGION 

Figure  1.  Electric  field  distribution  inside  an  IMPATT  diode. 
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1 


i „ m-l- 1 t , , ,,  m+1 

where  E,  lias  been  neglected  since  E.  << 
il  d 

field  at  breakdown  under  do  conditions  (i.e., 
stituting  Eq.  2 into  Eq . 1 yields 


m+ 1 
P 

when 


E is  the  peak 
po 

dx  = 1 ) . Sub- 


a dJ 
2 dt 


+ J 

s 


(3) 


The  peak  electric  1 ield  E can  be  written  as 

P 


v° 


E 

po 


E,  + E sin  nit 
b a 


t + t')  J (t')dt'  , (4) 


where  is  a constant,  E^sin  ait  is  the  field  from  the  applied  rf  voltage, 

T is  the  drift  zone  transit  time,  and  -1/a,  / (x,  - t + t')  J(t')dt' 

d d Jt-Td  d 

is  a field  resulting  from  the  so-called  space-charge  effect.  As  an 
approximation,  the  effects  of  E and  the  space  charge  field  can  be 
treated  by  assuming  that  the  carrier  generation  rate  is  sinusoidal  but 
has  a zero  at  tr  - 6 instead  of  at  tt.^  Therefore  Eq . 4 becomes 


E(t)  = E + E sin  (ait  + 5)  , 

po  a 

and  Eq . 3 can  be  reduced  to 


dd 

dt 


2 (m  + 1)  j _a_ 

T ' E 

a po 


sin(oit  + 6) 


(5) 


where  only  the  first  term 


of  the  expansion  of  (E  /E  ) 

P po 


m+1 


-1  was  kept. 


1 1 


B. 


FUNDAMENTAL  OPTICAL  INJECTION  LOCKING 


Suppose  the  IMPATT  diode  is  illuminated  by  an  optical  signal  modu- 
lated at  frequency  w;  then  .1  can  be  written  as 

.1  = .1  + J , sin  (cot  + 6 + (b)  . (6) 

s so  si 


Substituting  Eq.  b into  Eq . 5 yields 


dJ 

dt 


2 (m+l)E  2 

J — 3 sin  (ct  + 6)  + — 

fa  po  T a 


s i n ( t 


+ 6 + 


T 


Eq . 7 is  of  the  form 


(7) 


dy 

dt 


+ P(t)  y = Q(t ) 


(8) 


whose  solution  is 


y 


-Jp(t)dt  f J*P(t)dt  - J" P ( t ) d t 

e I Q (t ) e dt+Ce 


(9) 


where  C is  a constant. 

Comparing  Eqs.  7 and  8 gives 


2 (m  + 1)  E 

P(t)  = - — sin  (u)t  + 6) 

X E 
a po 

Q(t)  = — |J  + J , sin  (wt  + 6 + <|) ) 
' T I so  si 

a L 


and 


2 (m  + 1)  E 

JP(t)dt  = — — - cos  (wt  + 6)  = xcos  (wt  + 6) 


0)T  E 
a po 
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whero 


2 (m  + 1 ) 'a 


ui  E 

a po 


Therefore,  the  solution  for  J is  given  by 


-xeos (cot+b) 


jjJso  + Jsj  sin  (tot  + 6 + 4>)J  eXLOS  ^ ~ ~ ' dt 


; (tot+6) 


+ .1  e 
o 


-xcos (wt+6) 


(10) 


With  the  help  of  the  following  relations 


OO 

£ 


xcos(cot+6)  , . ^ _ 

e - I (x)  + 2 > I (x)  cosn(a)t  + 6) 

o / n 


n=] 


e-xcos(ait+(S)  = ^(x)  + (_i)n  j jx)  cos  n (tot  + 6)  , 

n=l 


n 


we  can  rewrite  Kq . 10  as 


J = J 


CX) 


I0(xH2  y ^(-1)  I ^ (x)cos  n(u)t+6) 


n=l 


oo 

£< 


(x)+2  In(x)cosn(ojt+6) 


n=  1 


Jc1 1 (x) 
si  o 


cos  (cot+6+(J>) 


j[JsoIo(x)  + Jsl1l  (x)sin(fl]  t - 

El  00  j . v-^  I (x)  r 

~n  sln  n(a,t+6)  - ~~  / y L cos  |^(n+l)  (tot+6  )+<|>J 


+ 


n=i 


n=l 


oo 

-^y^^TTT  cos  £(n~l ) (tot+6)  - *][  . 

n=2  ’ 


(11) 


1 
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By  expanding  the  summations  in  Eq . 10  and  keeping  terms  involving  modified 
Bessel  functions  of  order  0,  1,  and  2 only,  J becomes: 


I , (x)  I?(x) 

J = J,  - 2 J — — cos(o)t-HS)  + 2 J — 7 — — cos  2 (wt+6) 

dc  dc  I (x)  dc  1 (x) 

O O 

2 ( J I (s) 

+ — I (x)  <[J  I (x)  + J 1 (x)  sin4>]  t - — - — - cos(wt+6+40 

T o I so  o s 1 1 U) 

a I 

2 J I,  (x)  J I (x)  J I (x) 

H — sin(oJt+(>)  + - S-°  — sin  2(odt-H5)  - — — ' — cos[  2 (u)t+6  )+4>  ] 

j , f 7 (x ) J 1 7 ( x ) I 

- — — ^ — cos  [ 3(wt+6)+<}>]  4 cos  (wt+S-<f>)  I 

JCO  U)  1 

4 l,(x)  l J I (x) 

- {[J  I (x)4-J  , I . (x)  sin(j>]  t cos(wt+6) — — ° {cos  [2(wt+6)+4)]+cos<J>} 

T / so  o si  L Zu) 

a l 

J I (x)  .1  1 (x) 

■+ sin  2(u't+i!i)  4 [sin  3(wt+6)  + sin (wt4-6 ) ] 

j iMx)  J ,1  (x) 

- — — 7 — [ cosi  3 (wt+6)+<J>}+cos  (oJt+5+ij))] — [cost  4 (wt+6)+6} 

4U)  OCjO 

J ,I„(x)  ) 

+ cos  {2  (wt4-6)4-<J>}]  + — °-r [ cos{ 2 (wt+6 ) -$}  + cos<j>W 

2u  \ 

4 I (x)  i J I (x) 

+ <(J  I (x)+J  I (x)sim{)tcos2(cot+6)-  -s„  0 [cos{3(wt+6)+<{} 

T / so  o si  1 2gj 

a i 

J I (X)  J I (x) 

+ cos (wt+6-<}>)  ]4 [sin3(wt+6)-sir.(wt+6)  ] + — sin4(ojt+6) 

U)  ZCJ 

3 i.fx)  J I (x) 

- —7 [cos{4 (wt+6)+<J)}+cos(()] - — 7 [cos{5(a)t+6)+4>}+cos(u)t+6+4>)  ] 

4U)  b(i) 

j , I7 (x)  I 

+ - "2~ [ cos { 3 (wt+6 )+<})}  + cos  (wt+6+iJ ) ] | , (12) 


where  ,1 , J I (x)  is  the  dc  current  of  the  IMPATT  diode.  It  is  clear 
dc  o o 

from  Eq . 12  that,  although  the  applied  rf  voltage  across  the  IMPATT 
diode  is  sinusoidal,  the  corresponding  rf  current  is  not  sinusoidal. 
However,  we  are  interested  in  the  fundamental  component  only.  If  we 
denote  the  fundamental  component  of  J as  then 


J 


(1) 


I.(x)  J I (x)  [4  I (x)  - 6 I9(x)] 

-2  J - — cos  (ait  + 6)  + — — sin  (tot  + <S) 

dc  I (x)  cot 


+ r~2  I (x)  + I_  (x)  + — I0  (x)  cos  (cot  + 6 + <J>) 

m r I o 1 1 

a u 


3 2 


>]« 


Il(x)  , 


'dc  1 (:<)  )tOS 

o ' 


os  ( . t + ' ) - 


-’so  l0(x)[2I0(x)-3I2(x)]  _ 


dc 


0)T 


— sin(u)t+6)> 


+ 2 J 


si 


4 1 ^ (x ) - 61  2(x)  + 31,2(x) 
2 o 1 


cos  (uit  + 6 + <p)  . (13) 


Since  J .1  , , it  follows  that 

so  dc 


J 

so 


I (x) [21  lx)  - 3I,(x)] 

o o Z 


COT 

a 


« i 


We  can  rewrite  .1 


(1) 


as 


I,  (x) 

J.1S  = -2  J,  — — — r-  cos  (Mt  + 6 ) + 2 J 1 K(x)  cos  (uit  + 6 + (p)  , (14) 
(1)  dc  1 (x)  s si 


where 


. -1  Jso  '„<*>  - 3I2W1  . , 

0 = sin  — (-  o 

S J , CuT 

dc  a 


41-  (x)  + 31.  (x)  - 61  (x) 


K(x)  = 


6(ij 
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The  external  circuit  current,  which  is  obtained  by  propagating  J 
through  the  drift  region,  is  given  by 


(1) 


(15) 


To  estimate  the  locking  range,  a lumped  circuit  is  used,3  as  shown  in 

Figure  2.  Y , the  admittance  of  the  IMPATT  diode,  is  equal  to  G + iB 

11’ 

which  can  be  found  by  dividing  the  first  term  in  Eq . 15  by  the  rf  voltage 


V = -jV1  e 


and  adding  the  depletion  susceptance  ojC  : 

d 

w'[.  /un.  \ 

Y - . 2 sln  — JdcA  -n-r-O 

1 J wr  V,  I (x)  e /+  JwCd 


(16) 


The  load  admittance  is 


Y = (R  + jwL) 


-1 


2 2 2 

R + ui  L 


“ J 


toL 


2.2 


R + Ui  L' 


The  injection  current  IL  is  obtained  by  multiplying  the  second  term  in 
Eq.  15  by  the  area  of  the  diode: 

WTd 

sin  ~2~  j{u>(t-T  /2)+6+<J>} 

TL  = 2 K(X>  J-sl  A 6 • (17) 

d 

2 

where  A is  the  cross-sectional  area  of  the  diode. 
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The  circuit  equation  under  locking  is 


iL  = (Y,  + Yl)  V , 


(18) 


and  the  circuit  equation  under  the  free-running  condition  is 


W + W 


V . 


(19) 


In  the  following  analysis,  we  will  assume  that 


ail 


d , TT  c-  II 

T*K-  2 ' V*  2 


In  this  case,  Eq.  16  simplifies  to 


4 I1(X)  rdc 

i = " n TTxTvT  + j(i,cd  G1  + j B1 

o 1 


The  imaginary  part  of  Eq . 19  gives 


(0  1*  o o o t 

„ o . .2,22  L 

(i)  C,  - — r — — = 0 or  R + w L = — 
° d 2 2 2 o C, 

R + (i)  L d 


The  imaginary  part  of  Eq . 18  gives 


1 ( 1 T ) 

m L 


= coC 


toL 


d 2 2 2 

R + u>  L 


— (li)  + Ato'  C — — r- 

o d „2 


(to  + Ato)L 

o 


R + (w  + 2Ato  to  + Ato^)L^ 
o o 


(to  + Ato)  L 

“ (Wo  + Aw)  Cd  ‘ :2°'.  "2.2 
R + to  L 
o 


» 2Ato  C,  (to  C.L)  » 2Ato  C, 
d o d d 


1 - 


2Ato  to 
o 

2 2 2 
R + to  L 
o 


(20) 


(21) 


20 


where  it  is  assumed  that  the  series  inductance  resonates  with  the 
IMPATT  capacitance  at  the  operating  frequency  (i.e.,  u>2L  = 1).  From 
Eq.  17,  it  follows  that  I^Oj)  is  given  by 


1(1,)  = --  K(x)  1 sinij) 
mi,  n si 


Therefore, 


4K(x)  I&1 

2A,,)  = ■ ~nr~r~‘ sin<p 

i d 


The  locking  range  is  defined  as  the  maximum  2Au)  achievable,  and,  from 

Eq.  22,  (2Ato)  is  obtained  by  setting  sin0  = *1: 

M max 


Since,  from  Eq . 20, 


it  follows  that 


isl 

71  V1  Cd 


4 ¥*>  Tdc 

°1  IT  I (x)  V, 

o i 


IC(x)  I (x)"  wlr 
If(x) 


( 2 A<  o ) 


where 


= (wCdR) 


Eq.  24  shows  that  the  locking  range  is  proportional  to  the  optically 
generated  current  I and  is  inversely  proportional  to  the  circuit  Q 
and  the  bias  dc  current.  We  will  use  a few  numerical  examples  to 
estimate  (2Acj) 
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One  set  of  typical  values  is  presented  below: 


'dc 

= 25  mA 

E 

po 

2 x 107 

V/ m 

'si  " 

20 

PA 

l 

a 

= 1 

pm 

v = 

1 x 105 

m/ sec 

Q = 

20 

*d 

= 4 

Mm 

A = 

1 o"8  2 

1 0 m 

m = 

6 

= «a  + fd 

= 5 

Mm 

V1  = 

10  V 

f = 

10 

GHz 

Based  on  these  values,  it  follows  that: 

I = U /v  = 10  11  sec  an  =0.6) 

a a a 

Ea  = VL/fc  = 2 x 106  V/m 
2(m  + 1)  E 

x = — §.  = 2.22 

U)T  E z 

a po 

Io(x)  = 2.63  , 1 2 ( x ) = 1.91  , I2(x)  = 0.89 

K(x)  = -7.25 

(2Af)max  = 4,0  MHz  for  f = 1°  GHz. 

If  we  take  to  be  15  V rather  than  10  V,  then 

x = 3.33  , Io(x)  = 6.24  , ^(x)  = 5.23  , I2(x)  = 2.97, 

I<(x)  = -30.76 

and 


(2Af)max  = 14 • 7 MHz  for  f = 10  GHz. 


L 
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C.  FIRST  SUBHARMONIC  OPTICAL  INJECTION  LOCKING 
^WIMPATT:U>in  jection  2:1) 

Assume  that  the  optically  generated  current  is  at  the  first  sub- 
harmonic and  that  the  reverse  saturation  current  is  given  by 


= J + J . /n  sin  — (wt  +6+6) 
s so  sl/2  2 


Again  we  use  Eq . 5 to  solve  for  the  avalanche  current  J: 


J = J I 
o I o 


oo 

(x)  + 2 / (-1)'1  I (x)  cos  n (wt  +6)  + ~-  I (x) 

L-J  n Ta  ° 


*2U  -1)  In(x)  COS  n((0t  + 6)j  JsoIo(x)t 


2Jsl/2  Io(x)  I 


COS  7.-  (tot  + 6 + 4>) 


2 J 1 (x>  2J  (x) 

so  X n . , r.  sl/2  X n 

+ / — sm  n(<ot  + 6)  - / - — — r 

to  / j n to  / j 2n  + 1 


OO 

T.  ^TT005  [K)  (a)t+6)+f  1 


sl/2  \ * *n(x)  [7  1\ 

(0  2^J  2n  - 1 COS  li"  2/ 


(tot  + 6) 


Expanding  the  summations  up  to  n = 2 and  collecting  terms  involving 
cos  tot,  sin  tot , cos  1/2  tot,  and  sin  1/2  tot  yields 
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1 , (x) 

J = -2  J,  -1 


1 v ' 4 J / 

dc  TOO  cos  (<1,t  + «>  + — p ^ (U  j^(x)  ..  2 2 
o s wta  5 *2U;  3 Jj 

+ 2 - -3  <2<X>]  cos  ; 


(x)  - I 


H 


(wt  + 6 + rj,) 


+ 2 «i"4.  IjOofl^x)  -i,2(x,J  sin  i («  +6  + J 

_ I1(x) 

= ~2  Jdc  I~ 00  cos  (ut  +6  ) + j.  . _ 
o ' -s 


1/2 


wlie  re 


Jl/2  = A1  C°S  2 (a,t  + + *>  + A2  sin  I (Wt  + 6 + 0) 


A1  = 


4 Jsl/2  ( 


2 t2 


MT 


jl5  J2(x)  ~ 3 r^(x)  - To^>  + 2 cos ,p  I (X) 


a = Jsl/2  /->  r i 

2 ~ j2  Sln,/1  ri(x>  [io(x)  - j I2(x) 

and » therefore. 


fv»  - j >2m]I 


I 


Ji«-[*^]l/2>  Al 


/A2  + A 2 


cos  2 (ait  +6+0)  + 


A 


2 . .2]1/2 


]f[  + 


2 . 1 , , 
sinj  (ait  + 6 + 0)\ 


A1  + A2 


! sln  2 + <5  + <p  + 0) 


Jl/2o  sin  2 (wt  + 6 + <P  ) 


where 


<P  - 4>  + Q , and  tan  tt  = — 

2 A2 
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(26) 


n 


It  is  clear  that  0 is  a function  of  4>.  A plot  of  0 versus  <J>  is  shown  in 
Figure  3. 

The  avalanche  current  component  at  1/2  ait  (i.e.,  can  modify 
the  rf  electric  field  in  the  avalanche  region  through  the  space-charge 
effect.  The  modified  electric  field  is  E + AE  with  AE  given  by: 


AE(t) 


-±.L 


(Td  - t + t')  J1/7(t')  dt' 


Since  J^^.;(t')  varies  with  a period  of  4t(/o),  it  can  be  approximated  as  a 
constant  during  the  integration  time  interval: 


AE(t)  = - 2?  Ji/2(t) 


2T  J1/2o  sin  2 (ut  + 6 + V 


E,  sin  — (ait  + 6 + <\>  ) , 

D Z S 


where 


„ = Td  Jl/2o 

Lb  2t 


We  can  now  solve  for  the  avalanche  current  using  the  equation 


dJ  = 2 J (m  1) 
dt  T 


whose  solution  is 


sin  (ait  + 


po 


1 

— — sin  y (ait  + 6 + (Ji  ) 
*^?o 


J = J e 
o 


-xcos  (uit+6)+yc.osy  (wt+6-Kj)) 


(26) 


(27) 
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By  retaining  terms  involving  uit  only 


J = Jo  <Io(x)  I o ( y ) - 2 1 j (x)  cos  (oit  + cS) 


+ 2Iq(x)  l9(y)  cos  (tat  + A + t>^) 


+ 2 I (x)  1 2 ( y ) cos  (<i>t  + (S  - ,'i  j + "‘I 

( I , (x)  i2(y) 

.1  1 (x)  I (y)  <1-2  - — 7 — 7 cos  (cot  + 6)  + 2 -r—. — — cos  (oit  + 6 + (J-  ) 

o o o)  I (x)  I (y)  s 

\ o o 

l9(x)  I2(y)  | 

+ 2 I (x)  1 (y)  CO"  (,I,t  + S ‘ 'V  + '•  ( 

O O J 


= Jdc  - 2 Jdc  Too  cos  (wt  + :S)  + 2 Jdc  TTT 
o o 


cos  (ojt  + 6 + <p^) 


I9  (x) 

+ — — - — — COS  ((a)  t + iS  — (1)  ) 

Vx) 


The  external  circuit  current  is  then  given  by 


2 sin  — y~  I l^fx) 

J = J,  J , ~ cos 

e dc  on  , do  1 (x) 

d f o 


s [u^t  - -2d)  + <s] 


i.,(y)  ("  !2(x) 

J - — 1+2  — y cos  2 <p 

dc  I o ( y ) To(x)  s 


r2  n-iI/2 


l2(x)i  r / lA 

1 cos  “T  - Y) 

io(x)J  L v ’ 


+ 6 + cb  - ip 

s 


The  injection  current  is  given  by 


sin  ~y-  I (y)  T I 2 ( x ) 

I,  =21. T 1+2  7^7-r  cos  2 4> 

L dc  an  , 1 (y ) I (x)  s 

d ° L ° 


-il/2 

cos  ai^t + 6 + (J)^  - ip 


t 2 / .nl/2 
t .,  (x) 

+ — cos 

L (x). 
o 


If  we  assume  that  + 6 ~ ^ and  6 <<  ~ , then, 

from  F.q.  18,  we  can  write: 

*dc  4 *2^^  *2^x^ 

2 Aw  = cTX  v T77y  1 + 2 r~oo  cos  2 *s  + 77T  sin  (4,s  - ^ 

d 1 o 7 o I (x) 

^ o 


2 Aa)  = 


-1  2 ll/2 

I (x)  I,(y)  M I,(X)  I,(x) 

T / T't  7T  1+2  , -r  cos  2 cf>  + — r §in  (<J)  - iJO 

Tl(x)  Io(y)  Q IQ(x)  s X2(x) 

-*  o ' 


Therefore,  the  locking  range  is  the  maximum  value  of  the  right  side  of 
Eq.  29.  We  will  calculate  the  magnitude  of  2Aw  for  a number  of  p angles. 
The  following  numerical  parameters  are  used: 

I . = 20  pA  , 1 . = 25  mA  , V,  = 15  V . 

si/ 2 dc  l 

Q = 20  , = 1 pm  , II.  = 4 pm  . 

t =11  , E = 2 x 10^  V/m  , A = 10  ^ m^ 

o po 

v = 10^  m/sec  , f = 10  GHz 

x = 3.33  , lo(x)  = 6 . 24  , I j (x ) = 5.23  , l2(x)  = 2.97 

A^  = -7.11  x 10^  +6.97  x 10^  cos  <p 
A,(  = 6.97  x 10^  sin  (|i 

For 

(p  = 0°,  0 = 180°,  <J>  = -180°,  and  ip  = 0°  , 

we  have 

sin  (<f>s  - (|j)  = 0 and  2Aw  = 0 

For 

(p  = 10°  , 0 = -22.97°  , <p  = -12.97°  , ip  = -8.29° 
y = 0.06,  and  sin  = -8.16  x 10  ^ , 

we  have 

2Af  = 0.03  MHz. 
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Fo  r 


<J>  = 30 ° . o = -34.3°  , <}>s  = -4.3°  , ip  = -2.2° 

y = 0.166,  and  sin  (.[-^  - iJj)  = -0.04  , 

we  have 

2Af  = 0.12  MHz. 

Fo  r 

0 = 60°  , 0 = -61.97°  , j>  = -1.97°  , ip  = -1.27° 

y = 0.32  , and  sin  (<p  - ip)  = 1.27  x 10_2 

we  have 

2Af  = 0.124  MHz. 

For 

4,  = 90°  , (,  = -91.14°  , <Ps  = -1.14°  , <p  = -0.74° 

v = 0.45  , and  sin  (4>s  - i|»)  = 6.98  x 10"3  , 

we  have 

2Af  =0.12  MHz. 

from  the  above  numerical  values,  we  conclude  that  the  first  subharmonic 
optical  injection  locking  is  relatively  inefficient. 

D.  SECOND  SUBHARMONIC  OPTICAL  INJECTION  LOCKING 
(WIMPATT:a)injection  = 3:1) 

If  we  let  Jg  = sin  1/3  (tot  + 6 + <p)  and  solve  for  J 

in  Eq . 5,  then 
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E 

«=i 


J = -Jo  j I o ( x > + 2 P ( - 1 ) 11 1 n ( X ) cos  n(iut+6)j  + 1- 


(x) 


•X) 

E 


+ 2 7 (-1)  I (x)  cos  n( .it  + ) 

n=  1 


3.1. .,.1  (x) 

.1  l (x)t  - — t:os  (i.it  + <S  + 6) 

so  o J 


til 


E 

n=  1 


VI  (x) 

sin  » (<ut  + 6) 


'JJsl/3^  I '(3C) 


T.  3n  + 1 COS  (n  + i)(Wt  + ')  + 3] 

n=l 


+ ' Sl/>2^  in  - 1 COS  (n  - l)(  ' + ) ~ 3 


n=l 

Terms  involving  cos  1/3  (<ut  + 6 + <p)  are  grouped  to  give 
J 


si/  3 


1/3  cut 


a L 


12  2 3 2 2 

if  12(X)  ~ 2 Il(x)  " 6lo(x) 


cos  - (out  +6+41) 


= J^/3o  cos  "3  (aJt  + ^ + 4>) 


The  electric  field  is  modified  by  AE,  which  is  given  by 

AF  = ^ = F 

AE  2c  “ Ebl/3  • 
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Therefore,  the  new  r!  field  is 


- 2 1 ^ (y ) sin  (wt  + 6 + ij>)  + ••• 


= Jo  |‘o(x)  Vy)  ~ 2rc/y^  cos  ^'‘'t  + + 2Io(y)  I2(x)  cos  2 (o)t  + 6) 

+ 2Io(x)  I j (y)  sin  j (wt  + 6 + ip)  - 2Io(x)  I2(y)  cos  j (tot  + 6 + <J>) 

- 2Iq(x)  I^(y)  sin  (u>t  + 6 + <J>) 

-AI9(x)  I^(y)  cos  2 (wt  + 6)  sin  (tot  + 6 + <p)  + •••  > 
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It  wo  keep  only  terms  involving  .1,  then 


jd  v0(x)  ~ 2l(,(y)  * i(x)  °os  (“iL  + &) 


- 21^(x)  I }(y)  sin  (wt  + N + |) 


+ 2 1 , (x)  1 i ( V ) sin  (i.it  + 'S  - 'ji) 


>,(*)  i (y) 

J.  - 2.1  - --  cos  (u>t  + <S)  - 2.),  l 

<JC  do  I (x)  do  I (y) 

O ow 


/l2(x)' 
V I.  (*), 


2-il/2 


cos  (i.'t  + A + ip  - iji)  , 


2 — cos  2-(, 
o 


wlio  re 


tan  ip 


l2(x) 

1 ■ t~oo  cos  2 * 

o 

I2(x) 

Too  sin  2,|) 

o 


Therefore,  the  injection  current  is 


I (y)  sin  — i 

J,  =J1/_A=-2I.  ~ — <1 

h 1/3  dc  I (v)  art,  J 

o d 1 


U(x)  / 1 „ (x)\ 

2 17,0  cos  2*  + \Tw)  f 


2x1/2 


cos  1 0)  |t-  - + 4 + 4 - 4 1 


and  the  locking  range  is 


!o(x)  1 3 (y ) j ... 


^(x)  Io(y)J  Q 


I2(x)  I2(x)-)1/2 

2 fw  cos  2;  + ~2~~ 

O I (x) 


sin  (p  - 


14 


Use  the  following  numerical  values: 


Isl/3  = 

20  Vi  A , 

I . = 25  mA  , 

dc 

V1 

= 15 

V 

Q = 

20  , 

/ = 1 Jim, 

= 4 . 

.m 

a 

d 

_ 

li«  , , 

E “ 2 x 107  V/m  , 

A 

= 10' 

-8 

o 

po 

V = 

10  m/sec  , 

f = 10  GHz  . 

From  these,  we  get 


* = 3.33  , i()(x)  = 6.24  , ^(x)  = 5.23  , i2(x)  = 2.97 

-5  My) 

3]_/30  “ “8.62  x 10  , y - 0.6  , — — = 4.21  x 10 


For 


= 0°  and  'i  = 90°,  we  have  27. f = 1.32  MHz. 


For 


t-  = 15°  and  j = 68°,  we  have  2, If  = 1.27  MHz. 

For 

<P  - 30°  and  v = 61.6°,  we  have  2Af  = 0.99  KHz. 

For 


t-  = 60°  and  \ = 71.6°,  we  have  2.'.f  = 0.66  MHz. 
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DISCUSSION 


E. 

it  is  clear  from  the  above  calculations  that,  to  increase  the  locking 
range,  both  x and  y should  he  made  as  large  as  possible,  where 

2 (m  + I ) E 

a 

x ..H  ~ E 

a po 


2 (m  + 1 ) n J 


1/n  o d 


a po 


V , / . 


E is  the  rt  electric  field  in  the  diode.  Therefore,  in  a short  diode, 
a 

E^  is  larger  for  a given  . On  the  other  hand,  the  optimum  oscillation 
frequency  of  a shorter  diode  is  higher.  Thus,  the  optimal  situation 
seems  to  be  to  operate  the  diode  in  the  lower  portion  of  its  frequency 
band.  Our  limited  experimental  results  seem  to  indicate  that  this  is 
the  case.  The  best  way  to  increase  y is  to  have  large  optically  gen- 
erated current.  Numerical  values  used  in  the  previous  calculations  do 
not  represent  the  best  situation.  A simple  calculation  reveals  that 
1 mW  of  optical  power  at  rf  frequency  can  generate  up  to  660  uA  of  rf 
photocurrent  under  ideal  conditions.  This  would  mean  a locking  range  of 

250  MHz  for  the  case  of  w : u>  . . . = 3:1.  Therefore,  it  is 

osc  injection 

important  to  investigate  means  of  generating  efficient,  high-frequency, 
modulated  optical  sources. 

Preliminary  experiments  on  subharmonic  optical  injection  locking 
of  X-band  Si  IMPATT  oscillators  have  shown  a locking  range  of  close  to 
7 MHz  with  a 4 to  1 frequency  ratio.  The  photocurrent  (rf)  amplitude  in 
these  experiments  is  estimated  to  be  only  about  5 pA.  However,  the 
actual  numerical  parameters  of  the  diodes  used  are  not  available  to 
allow  comparison  with  analytical  results. 

Our  analysis  involves  several  approximations, and  the  result  is 
accurate  only  for  moderate  rf  voltage  swings  and  photocurrents.  Also, 
the  effect  of  circuit  tuning  is  completely  neglected.  We  believe  that 
in  a real  oscillator  the  interaction  of  the  device  and  the  circuit  can 


j 
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play  a major  role  in  determining  the  injection  locking  efficiency.  A 

more  rigorous  analytical  approach  is  being  taken  to  obtain  results  that 

are  applicable  under  even  large-signal  operations. 

It  is  interesting  that  injection  locking  with  an  input  signal 

proportional  to  sin  2/1  alt  is  also  possible.  J2/l  can  generate 

-xcoswt  ycos2/3u)t 

E,  “■  sin  2/1  alt,  and,  from  Eq . 27,  e e contains  a term 

b 

( r L (x)  cos  ait)  (I  j(y)  cos  2(ii  t)  that  produces  a current  component  at 
cos  ait . This  current  component  then  acts  as  the  driving  injection  sig- 
nal. Similarly,  any  multiples  of  the  subharmonics  are  possible  injection 
locking  signals  as  well. 

In  conclusion,  we  have  shown  that  subtiarmonic  optical  injection 
locking  of  Si  1MPATT  oscillators  can  have  a locking  range  of  0,100  MHz  if 
modulated  optical  signal  can  be  efficiently  converted  into  rf  photo- 
current. In  view  of  the  experimental  results  obtained  so  far  (more  than 
5 MHz  of  locking  range  with  less  than  5 pA  of  rf  optical  current),  we 
feel  that  the  calculated  results  underestimate  the  locking  range, 
although  the  experimental  device  parameters  are  not  known  well  enough 
to  allow  careful  comparison.  Nevertheless,  we  do  not  see  any  funda- 
mental reason  why  this  scheme  will  not  work  at  millimeter-wave 
f requenc ies . 
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SECTION  3 


RECENT  RESULTS  ON  OPTICAL  INJECTION  LOCKING 
OF  Si  IMPATT  OSCILLATORS 

Last  quarter  we  reported  some  initial  results  on  the  successful 
optical  injection  locking  of  X-band  Si  IMPATT  oscillators.  We  achieved 
locking  of  an  8.754-GHz  IMPATT  oscillator  by  injecting  an  optical  signal 
modulated  at  2.918  GHz.  We  also  observed  a reduction  of  oscillator  side- 
band noise  through  optical-illumination  and  injection-locking  processes. 

In  this  quarter,  we  continued  the  experiment  with  particular  emphasis  on 
obtaining  the  locking  range  versus  locking  gain  relation  for  various 
orders  of  subharmonic  injection  signal. 

The  injection- locking  characteristics  of  IMPATT  oscillators  were 
found  to  depend  heavily  on  the  cavity  tuning  conditions.  If  the  cavity 
was  tuned  such  that  the  IMPATT  oscillated  with  a high  Q factor  (i.e., 
narrow  spectral  width  and  minimum  side-band  noise),  the  locking  range 
typically  was  very  small  (a  few  hundred  kilohertz).  However,  if  the 
cavity  was  tuned  for  low-Q  oscillation,  then  the  locking  process  became 
very  efficient.  This  observation  agreed  well  with  our  theoretical 
calculations,  given  in  Section  2.  Another  factor  that  affected  the 
injection  locking  process  was  the  dc  bias  of  the  IMPATT  diode.  We  found 
that,  depending  on  the  bias  condition  of  the  diode,  different  amounts  of 
IMPATT  oscillation  frequency  shift  could  be  observed  for  a fixed  dc 
optical  illumination.  Under  certain  conditions,  the  frequency  shift  was 
positive;  but,  for  different  combinations  of  circuit  tuning  and  dc  bias, 
it  could  also  be  negative.  In  general,  a large  locking  range  was  obtained 
if  the  oscillator  frequency  shift  due  to  optical  illumination  was  large. 

One  of  our  IMPATT  oscillators,  when  biased  at  28  mA  and  tuned  to 
8.116  GHz,  could  be  optically  injection  locked  to  several  different  orders 
of  subharmonics.  Figure  4 shows  the  changes  of  oscillator  spectrum  at 
five  different  injection  signal  power  levels.  From  Figure  4(a)  to 
Figure  4(d),  the  injection  signal  level  was  raised  from  -20  dBm  to  10  dBm 
in  10  dBm  steps.  As  shown  in  Figure  4(a),  the  oscillator  output  was 
initially  at  8.111  GHz  while  the  injection  signal  was  at  roughly 
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Figure  4 


Spectra  of  an  injection-locked  [MPA IT  oscil 
lator  at  various  injection  signal  levels. 
1MPAYT  oscillation  frequency:  8.111  OHz 
Injection  signal  frequency:  1 . 152  (ill/ 
Horizontal  scale:  2 MHz/div 
Vertical  scale:  10  tlB/div 
Reference  level : JO  dBm. 


one-sixth  of  8.1146  GHz,  or  1.3524  (.'l!x.  However,  as  shown  in  Figure  4(e), 
when  injection  locking  finally  took  place,  the  oscillator  frequency  was 
pulled  in  and  locked  to  the  sixth  harmonic  of  the  injected  signal 
f requency . 

Figure  5 is  a plot  of  the  locking  band  versus  locking  gain  of  the 
same  IMPATT  oscillator  for  several  different  subharmonic  orders.  The 
1MPATT  oscillation  frequency  was  at  8.116  GHz.  A frequency  ratio  as 
high  as  8:1  was  used  (i.e.,  the  injection  signal  was  running  at  1.014  GHz) 
to  achieve  injection  locking.  The  locking  gain  was  defined  as  the  ratio 
of  the  IMPATT  oscillator  output  power  to  the  microwave  power  used  to 
modulate  the  injection  laser.  Therefore,  the  plot  of  locking  band  versus 
locking  gain  actually  contains  the  information  on  modulation  response  of 
the  injection  laser  as  well.  For  instance,  the  modulation  depth  of  the 
laser  at  2.03  GHz  (4:1  frequency  ratio)  was  not  as  efficient  as  at 
1.35  GHz  (6:1  frequency  ratio);  therefore,  the  injection  locking  efficiency 
of  the  6:1  subharmonic  was  better  than  that  of  the  4:1  subharmonic  when 
the  locking  gain  was  greater  than  8 dB.  It  is  also  interesting  that  in 
our  oscillator  the  even  order  of  subharmonics  is  preferred  to  the  odd 
order  of  subharmonics  with  respect  to  the  injection  locking  efficiency. 

We  believe  that  Lhis  is  a resulL  of  the  particular  cavity  tuning. 

Further  studies  are  needed  to  clarify  this  point. 
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Figure  5.  Locking  band  versus  locking  gain  of  a subharmonically 
optically  injection-locked  IMPATT  oscillator. 
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SECTION  4 


MODE  LOCK  I NO  OK  INJECTION  LASERS 


In  the  last  report,  we  discussed  the  genertil  effects  of  laser 

medium  dispersion,  laser  relaxation  oscillation,  and  the  laser  spectral 

broadening  mechanism  on  the  success  of  mode-locking  GaAs  injection 

lasers.  The  first  demonstration  of  injection  laser  mode  locking  at 
4 

MIT  was  achieved  by  using  an  external-cavity  cont iguration  and  by  vary- 
ing the  current  through  the  laser  diode  at  a rate  determined  by  the 
photon  round-trip  transit  time  in  the  cavity.  This  section  expands  on 
that  discussion  with  emphasis  on  the  effects  of  different  external 
cavity  configurations. 

One  of  the  reasons  for  using  an  external  cavity  for  injection  laser 
mode  locking  is  to  reduce  the  required  gain  modulation  frequency  to  a 
few  gigahertz.  Without  an  external  cavity,  the  required  frequency  f is 
80  to  130  GHz  for  typical  diode  lengths  (300  to  500  pm),  where 


c is  the  speed  of  light,  n is  the  index  of  refraction  of  GaAs,  and  i.  is 
the  diode  length.  Since  the  modulation  response  of  the  laser  decreases 
sharply  beyond  a few  gigahertz,  it  would  be  extremely  difficult  in  the 
80-  to  130-GHz  range.  Theoretically,  it  is  possible  to  achieve  laser 
mode-locking  witli  infinitesimal  perturbation  if  the  driving  frequency  is 
exactly  equal  to  the  transit  time  frequency.'’  Practically,  however, 
there  are  several  limitations.  First,  there  will  always  be  some  detuning 
caused  by  mechanical  vibrations  and  temperature  fluctuations.  Second, 
the  phase  noise  in  the  signal  source  and  in  the  laser  output  will  also 
effectively  detune  the  driving  frequency.  Thus,  there  is  a minimum 
perturbation  needed  for  mode  locking.  This  sets  a maximum  frequency 
usable  for  mode  locking  an  injection  laser.  By  using  an  external  cavity, 
the  problems  associated  with  very  high  frequency  laser  modulation  can  be 
ci  rctimvented . 
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Another  reason  for  using  the  external  cavity  configuration  is  that 
the  effects  of  material  dispersion  can  be  reduced  signil  iiantly.  This 
is  because  ol  the  Large  dispersion-free  air  space  in  the  optical  cavity, 
which  renders  the  gain  medium  dispersion  insignil icat.t  since  it  con- 
stitutes only  about  1%  of  the  total  cavity  length.  The  material  dis- 
persion effect  is  important  in  an  injection  laser  without  external 
cavity  since,  in  that  case,  the  longitudinal  modes  ot  the  laser  are 
spaced  unequally.  However,  under  mode-locked  conditions,  the  mode 
spacings  of  the  locked  longitudinal  modes  are  equal.  That  means 
that  the  gain  perturbalion  should  be  strong  enough  to  pull  the  modes  to 
their  proper  frequency  Locations  in  order  to  lock  them  together.  Again, 
this  will  set  a limit  on  the  minimum  input  driving  power  required  to 
obtain  a given  mode-locked  laser  pulsewidth  (which  is  inversely  propor- 
tional to  the  number  of  modes  locked)  or  on  the  minimum  pulsewidth 
obtainable  for  a given  input  driving  power. 

To  induce  coupling  between  longitudinal  modes  of  a laser  requires 
a gain  (or  loss)  modulation  with  spatial  dependence.  The  coupling  coef- 
ficient between  two  longitudinal  modes  can  be  expressed  as^ 

ka , b = /jJvA8(x,y,Z^  Ea(x»y>z)  F.b(x,y,z)  dxdydz  , 

where  E and  E,  are  two  normal  modes  of  the  laser  cavity,  Ag  is  the  gain 
a d 

variation,  and  the  integral  is  over  the  cavity  volume  v.  Since  E and 

3 

E^  are  orthogonal,  i.e., 

Ufy  E a • Eb  dxdydz  = 0 ’ 

Ag  should  be  a nonuniform  perturbation  if  the  coefficient  of  mode 
coupling  is  to  be  nonzero.  The  use  of  an  external  cavity  provides  a 
convenient  solution  for  this  problem. 
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There  are  two  external-cavity  configurations  that  can  be  used  for 
the  injection  laser  mode-locking  experiment.  We  will  examine  their 
transmission  characteristics  and  determine  which  configuration  is  more 
suitable  for  our  purpose.  The  two  configurations  are  shown  in  Figure  6. 
In  Figure  6(a),  the  cavity  is  formed  by  one  cleaved  iacet  (A)  of  the 
semiconductor  laser  and  an  external  mirror  with  high  reflectivity 
(r  = 0.98).  An  AR  coating  is  applied  to  the  other  facet  (B)  of  the 
semiconductor  laser.  Figure  6(b)  shows  a more  complicated  external 
cavity,  which  is  essentially  the  same  as  the  first  cavity  except  that 
there  is  no  AR  coating  on  facet  B.  The  transmission  characteristics  of 
these  resonators  can  be  analyzed  by  regarding  the  composite  resonators 
as  one  simple  Fabry-Perot  resonator  of  length  1 . One  of  the  mirrors  of 
this  equivalent  resonator  is  mirror  C,  the  other  is  the  semiconductor 
laser.  The  semiconductor  laser  itself  is  a Fabry-Perot  of  length  £ ; 
its  equivalent  transmission  and  reflection  coefficient  can  be  calculated 
easily.  Thus,  the  composite  resonator  transmission  coefficient  is 
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Figure  6. 


Two  different  external-cavity  arrangements  for 
semiconductor  lasers. 


Eq . 30  gives  the  transmission  characteristics  of  the  cavity  shown 

in  Figure  6(b).  If  we  use  parameters  typical  of  the  GaAs  laser  and  let 

r j = 0.98,  Eq . 30  can  be  plotted  as  shown  in  Figure  7.  It  exhibits  the 

fine  longitudinal  mode  structure  of  a cavity  with  length  + 1‘  ~ 0, 

o I o 

and  the  beating  effect  of  the  two  cavities  of  length  y.  and  Si  + S’  . 

o o 1 

The  beat  frequency  is  c/2ii'y  If  we  set  r = 0 and  T^  = 1 in  Eq.  30, 
we  can  also  plot  the  transmission  characteristic  of  the  cavity  shown  in 
Figure  6(a).  This  is  displayed  in  Figure  7.  Since  the  cavity  structure 
is  basically  a simple  Fabry-Perot,  the  familiar  equally  spaced,  equal 
height  transmission  pattern  is  obtained,  as  shown  in  Figure  8.  How- 
ever, the  cavity  Q is  much  lower  in  the  simple  structure.  This  can  be 
seen  by  comparing  the  height  and  width  of  the  transmission  peaks  given 
in  Figures  9 and  10.  We  believe  that  the  three-mirror  cavity  is  a 
better  configuration  for  the  mode-locking  process  because  of  the  higher 
Q factor  and  the  additiomil  beating  effect  (which  forces  longitudinal 
modes  to  be  locked  in  groups).  Therefore,  the  locking  process  would  be 
more  stable  in  this  case. 

Experimental  work  is  underway.  An  external  cavity  injection  laser 
has  been  set  up.  We  have  observed  laser  threshold  current  reduction  due 
to  the  addition  of  the  external  mirror.  More  experimental  results  will 
be  presented  in  the  next  quarterly  report. 
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Figure  7 


Intensity  transmission  characteristics  of  the 
three-mirror  cavity  depicted  in  Figure  6(b). 
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Figure  9. 


Resonance  transmission  peak  of  the 
simple  cavity  shown  in  Figure  6(a). 
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Figure  10.  Resonance  transmission  peak  of  the  three-mirror 
cavity  shown  in  Figure  6(b). 
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SECTION  5 


PLANS  FOR  THE  NEXT  QUARTER 

In  the  next  quarter,  we  will  continue  the  injection  laser  mode- 
locking experiment.  Our  initial  results  indicate  that  the  external 
cavity  could  lower  the  laser  threshold,  although  the  alignment  would  be 
very  critical.  We  will  investigate  a slightly  different  arrangement. 
Instead  of  using  a highly  reflective  external  mirror,  we  will  use  a 
semi-transparent  mirror  such  that  the  transmitted  optical  power  is  the 
output  of  our  resonator.  Since  this  approach  will  eliminate  the  need 
for  specially  mounted  laser  diodes, the  experiment  can  be  carried  out 
with  lasers  that  have  only  one  emitting  facet  accessible.  Several 
optical-microwave  devices  will  be  analyzed  and  fabricated.  The  fabri- 
cation of  waveguide-incorporated  GaAs  FETs  is  underway.  The  devices 
will  be  tested  for  their  improvement  in  optical  response.  Major  efforts 
will  be  focused  on  the  study  of  high-speed  optical-detection  techniques. 
Our  approach  is  to  take  advantage  of  the  speed  of  microwave  solid-state 
devices  and  the  bandwidth  (and  possibly  the  gain)  derived  from  the 
oscillator  and  amplifier  circuits. 
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